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We present a systematic study of the B meson light-cone wavefunctions in QCD in the heavy-quark limit. We 
construct model-independent formulae for the light-cone wavefunctions in terms of independent dynamical degrees 
of freedom, which exactly satisfy the QCD equations of motion and constraints from heavy-quark symmetry. The 
results demonstrate novel behaviors of longitudinal as well as transverse momentum distribution in the B mesons. 


Recently systematic methods based on the 
QCD factorization have been developed for the 
exclusive B meson decays into light mesons [ Q 
(see also Ref. [|^). Essential ingredients in this 
approach are the light-cone distribution ampli¬ 
tudes for the participating mesons, which express 
nonperturbative long-distance contribution to the 
factorized amplitudes. The light-cone distribu¬ 
tion amplitudes describe the probability ampli¬ 
tudes to find the meson in a state with the con¬ 
stituents carrying definite light-cone momentum 
fraction, and thus are process-independent quan¬ 
tity. For the light mesons (tt, K, p, K*, etc.) 
appearing in the final state, systematic model- 
independent study of the light-cone distributions 
exists for both leading and higher twists [ §| . On 
the other hand, the light-cone distribution am¬ 
plitudes for the B mesons are not well-known at 
present and they provide a major source of un¬ 
certainty in the calculations of the decay rates. 

By definition, the distribution amplitudes are 
obtained from light-cone wavefunctions at (al¬ 
most) zero transverse separation of the con- 
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stituents, ^ ^^2 d^kT^{x,kT)- The 

light-cone wavefunctions with transverse momen¬ 
tum dependence are also necessary for computing 
the power corrections to the exclusive amplitudes, 
and for estimating the transition form factors for 
B ^ D, B ^ TT, etc, which constitute another 
type of long-distance contributions appearing in 
the factorization approaches for the exclusive B 
meson decays. 

In this work [ |^, ||, we demonstrate that, in 
the heavy-quark limit relevant for the factoriza¬ 
tion approaches for the exclusive B meson decays, 
the B meson light-cone wavefunctions obey exact 
differential equations, which are based on heavy- 
quark symmetry and the QCD equations of mo¬ 
tion. As solution of those differential equations, 
we derive the model-independent formulae for the 
light-cone wavefunctions, which involve not only 
the leading Fock-states with a minimal number of 
(valence) partons but also the higher Fock-states 
with additional dynamical gluons. 

The light-cone wavefunctions are related to 
the usual Bethe-Salpeter wavefunctions at equal 
light-cone time z~^ = + z ^)/In the heavy- 

quark limit, the quark-antiquark light-cone wave- 
functions ‘&±(t, z^) of the B mesons can be intro¬ 
duced in terms of vacuum-to-meson matrix ele¬ 
ment of nonlocal operators in the the heavy-quark 
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effective theory (HQET) [ 0 , |0: 


{0\q{z)Thy{0)\B{p)) = 2 Tr 

-pl+ 
75r 2 




■( 1 ) 


Here = (0, z~,zt), = —zip, = l,t = vz, 

and = Mv^ is the 4-inomentum of the B me¬ 
son with mass M. hy{x) denotes the effective 
6-quark field, b(x) « exp(—• x)hy{x), and 
is subject to the on-shell constraint, = hy 
[ 1^. r is a generic Dirac matrix and, here 
and in the following, the path-ordered gauge 
factors are implied in between the constituent 
fields, fs is the decay constant defined as usually 
as ( 0 |g( 0 ) 7 '' 75 /i„( 0 )|i?(p)) = ifsMv^, so that 
= 0, = 0) = 1. Eq. (P is the most gen¬ 

eral parameterization compatible with Lorentz in¬ 
variance and the heavy-quark limit. 

Higher Fock components in the B mesons are 
described by multi-particle wavefunctions. We 
explicitly deal with quark-antiquark-gluon three- 
particle wavefunctions, defined as [ |j 


{0\q{z) gGfj,,,{uz) z^T hy{0)\B{p)) 


75r 

u) - 4'y (t, - i a^yz'' u) 


= \fBMTr 


- Zfj,XA{t,u) + ^YA{t,u) 


.( 2 ) 


where the ellipses stand for the terms which in¬ 
volve one or more powers of z^ and are irrelevant 
for the present work. We have the four functions 
^v,^a,Xa and Ya as the independent three- 
particle wavefunctions in the heavy-quark limit. 

The QCD equations of motion impose a set of 
relations between the above wavefunctions [ii- 
They can be derived most directly from the exact 
identities between the nonlocal operators: 


—q{x)-t^Thy{Q) 


= I 


dun q{x)gG^^rhy{0) , (3) 




= i du{u — 1) q{x)gGfj,jp{ux)v^x’^rhy{0) 

Jo 

( 4 ) 

y^o 

where G/^y is the gluon held strength tensor, and 
we have used the equations of motion pq = 0 and 
V ■ Dhy = 0 with — igAf^ the covariant 

derivative. Taking the matrix element with ^ 
z^, the LHS of these identities yield l>+(t, z^), 
$_(t, z^) dehned in Eq. (j^) and their derivatives, 
d^±{t, z‘^)/dt and z^)/9z^. The terms in 

the RHS, which are given by integral of quark- 
antiquark-gluon operator, are expressed by the 
three-particle wavefunctions of Eq. (||) . The last 
term of Eq. (^), the derivative over the total 
translation, yields $±(6, z^) multiplied by 




A = M — mb = 


iv ■ d{0\qThy\B{p)) 
{0\qrhy\B{p)) 


( 5 ) 


This is the usual “effective mass” of meson states 
in the HQET [0. 

Substituting all the Dirac matrices for T, we 
obtain the four independent constraint equations 
between the two- and three-particle wavefunc¬ 
tions from Eqs. (i and (!)■ We solve this sys¬ 
tem of equations for the relevant two cases: (i) In 
the light-cone limit z^ —> 0 {zt —> 0), with fully 
taking into account the contribution due to the 
three-particle wavefunctions. The solution gives 
exact model-independent representaions for the 
light-cone distribution amplitudes in terms of in¬ 
dependent dynamical degrees of freedom, (ii) For 
z^ P 0 but in the approximation neglecting the 
contribution of the three-particle wavefunctions. 
The solution gives exact analytic formulae for the 
light-cone wavefunctions with transverse momen¬ 
tum dependence within the valence Fock-states. 

Now we discuss the case (i) in detail [0. In 
the light-cone limit, the light-cone wavefunctions 
of Eq. (0) reduce to the light-cone distribution 
amplitudes (t) as [ 0 

^±it) = = 0) , ( 6 ) 


and we also introduce the shorthand nota¬ 
tions, = d^±{t)/dt and d^±{t)/dz'^ = 

z^)/9z^|22^0j which denote the deriva- 
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tives with respect to the longitudinal and trans¬ 
verse separations, respectively. The first identity 
(||) yields the two equations: 


4)- 


= 2t duu(^A{t,y') — ^v{t,y')), (7) 

^0 

6+(t) - - I ((j)+{t)- 


t 


dz^ 


= 2t duuA A{t,u) + 2^vit,u) -I-u)), (8) 
Jo 


and similarly the second identity yields 

^+(0 - ^ (<^+(0 - ^-(O) + 


dz^ 


= t du{u-l){'SA{t,u) + XA{t,u)) , A) 

Jo 


b'+iO - J-iO + ( j [A{t) - <t>- 


+2t 


( dA{t) _ d^-jt) 

5z2 dz^ 


= 2t du {u — 1) AA{t,u)+ YA{t,u)).{W) 
Jo 

These Eqs. (0)-(^ are exact in QCD in the 
heavy-quark limit. 

Important observation is that we can eliminate 
the term dA{t)/dz'^ by combining Eqs. and 
(3). The resulting equation, combined with Eq. 
(7), gives a system of two differential equations 
which involve the degrees of freedom along the 
light-cone only. By going over to the momen¬ 
tum space by (j)±{t) = J duj e~'^‘^*(j)±{uj), and 
F{t,u) = J dujd^ with F = 

{dty, Xa, TaI, the corresponding differential 
equations read [ Q 

( 11 ) 

{uj — 2K) + ui(j)-{ijj) = J(ijj) , (12) 

where /(w) and J(tt’) denote the “source” terms 
due to three-particle wavefunctions as 


d /■“ 


dC d 


J{uj) = -2 


-4 


dio 


dp 


dp 


d^d^v{p,0 


[^A{p,0+XA{p,i)] 

(13) 


.-p C di ■ 

Eqs. (0). (0) can be solved for (/>+(a;) and 
4>-{u)). The boundary conditions are specified 
as 4>±{io) = 0 for ui < 0 or tu ^ oo, be¬ 
cause ujv'^ has the meaning of the light-cone pro¬ 
jection of the light-antiquark momentum in 
the B meson, and the normalization condition is 
du}(j)±{oj) = $±(0,0) = 1. Obviously, the so¬ 
lution can be decomposed into two pieces as 

4>±{u}) = (w) + (w) , (14) 

where (uj) are the solution with Hoj) = 
J{uj) = 0, which corresponds to the “Wandzura- 
Wilczek approximation Ty = 'I'a = Xa = 
Ya = 0. denote the pieces induced by the 

source terms I{uj) and J(w). 

We are able to obtain the analytic solution for 
the Wandzura-Wilczek part as 




A i (a; — A) 
2A2 


6»(w)6»(2A - w) . (15) 


Moreover, the solution for can be obtained 
straightforwardly, and reads (w > 0): 


^+\^) = ^^(^) 


b^f\u;) = 


2A — Lu 


GiA 


J{uj) 


2A 

g{uj) = 0{2A-uj) 
-0{uj - 2A) 


dp 


dp 


Kjp) 
2A- p 
K{p) 


(16) 

(17) 

2A I 


2A-P 

- 

with K{p) = I{p) + [l/(2A) — d/dp\ J{p). The 
solution dl^ ) with Eqs. 0)-® is exact, and 
reveals that (j)± contain the three-particle contri¬ 
butions. This is also visualized explicitly in terms 
of the Mellin moments (w"')± = f doj u;"0±(±;) 
(n = 0,1, 2, • • •). Here we present some examples 
for a few low moments: (uj}± = (3 ± l)A/3, and 


2\ 4±2^2 

n± = A 


1 \2 
3^^ 


^ \ 2 




E 5 


(19) 
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where Xe and Xh are due to and are 

related to the chromoelectric and chromomag- 
netic fields in the B meson rest frame as 
{0\qgE ■_a-f 5 hy\B{p = 0)) = fBMX%, {0\qgH ■ 
cr'-f 5 hv\B{p = 0)) = ifBMXfj. Our solution ([l^ ) 
allows us to further derive the analytic formulae 
for the general moment n in terms of matrix el¬ 
ement of local two- and three-particle operators 
with dimension n-|-3 (see Ref. for the detail). 

The behavior of Eq. for a fast-moving 

meson, t = v ■ z ^ oo, shows that </)+ is of 
leading-twist whereas (p- has subleading twist; 
the three-particle contributions to the leading- 
twist (pE s-re in contrast with the case of the 
light mesons [ where the leading-twist ampli¬ 
tudes correspond to the valence Fock component, 
while the higher-twist amplitudes involve multi¬ 
particle states. We note that there exists an esti¬ 
mate A|/A2 = 0.36 ± 0.20, A|^/A2 = 0.60 ± 0.23 
by QCD sum rules [ ^ (see Eq. (H))- This 
might suggest that, in the B mesons, the three- 
particle contributions could play important roles 
even in the leading twist level. In this connec¬ 
tion, the shape of our Wandzura-Wilczek con¬ 
tributions (|l5|), which are determined uniquely 
in analytic form in terms of A, is rather dif¬ 
ferent from various “model” distribution ampli¬ 
tudes that have been used in the existing litera¬ 
ture: One exapmle of such models is {to) = 
and {lj) = (l/wo)e-“/“'> with 
u)Q = 2A/3, inspired by the QCD sum rule es¬ 
timates [ These have very different shape 
compared with Eq. except the behavior 

~ w, (jP^{Lu) ~ const, as w ^ 0. 

Next we proceed to the case (ii), where we get 


o d 

a;^ + <I>+ + z2^(ci>+-ci>_) = 0, (20) 


duj 

d 


dz^ 


^^ + 2)(‘i>+-‘i>-)+4—^=0, (21) 


/ T, 9 3' 


1 

$ — 


92 3$+ 


2 ' 9^2 

($+-$_) 


= 0 ,( 22 ) 


+ 2 - 


/ 9$ 


a4> 


= 0 , (23) 


9^2 V az2 

corresponding to Eqs. (g), (|), (|), (^, 


re- 


spectively. Here Eqs. (|^-(|3|) are given in 
the “^-representation” = $±(0;, z^)^ in¬ 

stead of the “t-representation”, via $±(t, z^) = 
f doj e“®“*$±(a;, z2). The light-cone limit is not 
taken so that the terms proportional to ap¬ 
pear in Eq. (p^. Note that we have neglected 
the contribution from the quark-antiquark-gluon 
three-particle operators. 

By combining Eqs. (ED and d^), we eliminate 
the last term in their LHS. The resulting equa¬ 
tion can be integrated with boundary conditions 
$±(0;, z2) = 0 for tu < 0 or tu —> 00 as 

(w - 2A) 4>+ -f = 0 . (24) 

In the limit z^ ^ 0, Eqs. ( ^ and (|^ reduce 
to Eqs. and dll) in the Wandzura-Wilczek 
approximation. The corresponding solution 
serves as “boundary conditions” to solve Eqs. 
(2C)-(|^) for z2 7^ 0. Then, from Eqs. (|^ ) and 
(24), we find, as the solution in the Wandzura- 
Wilczek approximation for z^ 7^ 0, 

^^^\ijj,z^) = (j)^^\Lo)^[z'^ui{2K —uj)) . (25) 

Here ^(a;) is a function of a single variable 
X, and can be determined from a remaining 
differential equation (^ or ( p^ ) as ^(x) = 

Jo {V~^) where Jq is a (regular) Bessel function. 
This result gives analytic solution for the light- 
cone wavefunctions with the transverse separa¬ 
tion Zy = — z^. The momentum-space wavefunc¬ 
tions <i)£^^(aj, fcr), defined by —z|.) = 

/ dtod'^kx fey), read 

AW), X 

kr) = — - S (fe^ — a;(2A — tu)) .(26) 

The result (|^ gives exact description of the va¬ 
lence Fock components of the B meson wave- 
functions in the heavy-quark limit, and represent 
their transverse momentum dependence explic¬ 
itly. These results show that the dynamics within 
the two-particle Fock states is determined solely 
in terms of a single nonperturbative parameter A. 

The transverse momentum distributions in the 
B mesons have been unknown, so that various 
models have been used in the literature. Fre¬ 
quently used models assume complete separa¬ 
tion (factorization) between the longitudinal and 
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transverse momentum-dependence in the wave- 
functions (see e.g. Refs. [ II)- A typical exam¬ 
ple of such models [||] is given by fey) = 

^ x ) with some 

constants N, (3, and k. Eq. ( p6| ) shows that trans¬ 
verse and longitudinal momenta are strongly cor¬ 
related through the combination fc^/[a;(2A — w)], 
therefore the “factorization models” are not jus¬ 
tified. We further note that many models assume 
Gaussian distribution for the fer-dependence as 
in . These models show strong dump¬ 

ing at large \zt\ as ~ exp (^—k'^z'^/2). In con¬ 
trast, our wavefunctions (|^ ) have slow-damping 
with oscillatory behavior as —2^) ~ 

cos(|2:T|\/a;(2A - w) - 7r/4)/y^|^. 

Finally, we can estimate the effects neglected 
in our solutio n (|^ ). Inspecting the t —> 0 limit 
of Eqs. (^), (|1C|), one immediately obtains the 
exact result for d(j)±{0)/dz‘^, which gives the first 
moment of as [ ^ 

Jdud^kT 4$±(a;,feT) = I (A^ + A| + A|,), (27) 


with Xe and Xh of Eq. (|l^). Here -|- 

denote the total wavefunctions which include 
the higher Fock contributions induced by the 
three-particle operators. From Eq. (^), we get 
J duxPkT kr) = 2A^/3, so that the 

terms 2(A|; -I- A|^)/3 of Eq. (|^ come from 
The result (^), combined with a QCD sum rule 
estimate of Xe,Xh mentioned below Eq. (^, 
suggests that the higher Fock contributions might 
considerably broaden the transverse momentum 
distribution. However, qualitative features dis¬ 
cussed above, like non-factorization of longitudi¬ 
nal and transverse directions, “slow-damping” for 
transverse directions, etc., will be unaltered by 
the effects of multi-particle states. 

To summarize, we have derived a system of 
differential equations for the B meson light-cone 
wavefunctions and obtained the corresponding 
analytic solution. The differential equations are 
derived from the exact equations of motion of 
QCD in the heavy-quark limit. The heavy-quark 
symmetry plays an essential role by reducing the 
number of independent wavefunctions drastically, 


so that the configuration of quark and antiquark 
in the B mesons is described by only two light- 
cone wavefunctions. As a result, a system of 
four differential equations from the equations of 
motion allows us to obtain model-independent 
formulae of these two wavefunctions, which re¬ 
veal roles of the leading Fock-states, as well as 
the higher Fock-states with additional dynam¬ 
ical gluons. Also due to the power of heavy- 
quark symmetry, our Wandzura-Wilczek parts 
dl^), (p^), which correspond to the leading Fock- 
states, are given in simple analytic formulae in¬ 
volving one single nonperturbative parameter A. 
Heavy-quark symmetry also guarantees that our 
solutions determine the light-cone wavefunctions 
for the B* mesons and also for the D, D* mesons 
in the heavy-quark limit. 

We emphasize that our solutions provide the 
powerful framework for building up the B meson 
light-cone wavefunctions and their phenomeno¬ 
logical applications, because the solutions satisfy 
all relevant QCD constraints. Further develop¬ 
ments like those required to clarify the effects of 
multi-particle states can be exploited systemati¬ 
cally starting from the exact results in this work. 
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